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Recently, a method for 0(a) improvement of composite operators has been proposed wh ich uses the larg e 
momentum behavior of fixed gauge quark and gluon correlation functions (G. Martinelli el ai, bep-lat/0106003). 
A practical problem with this method is that a particular improvement coefficient, cngi, which has a gauge non- 
covariant form, is difficult to determine. Here I work out the size of the errors made in improvement coefficients 
and physical quantities if one does not include the cngi term. 



Non-perturbative renormalization jlj] (NPR) is 
a method for determining the renormalization 
constants for composite operators which does not 
rely on perturbation theory. In a recent paper, 
the method was extended to allow the determi- 
nation of the improvement coefficients needed to 
remove all 0(a) corrections from matrix elements 
of bilinear operators This improvement to 
NPR is important because it allows the determi- 
nation of fully 0(a) improved and renormalized 
composite operators even for operators with non- 
vanishing anomalous dimensions, e.g. scalar and 
pseudoscalar bilinears. 

The NPR method uses amputated quark and 
gluon correlation functions, determined in a fixed 
gauge, usually Landau gauge. To improve the 
method one must first improve the quark field 
itself, which, it is argued in Ref. H, requircsn 
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Here W =fl+mo is the entire O(a) improved Wil- 
son action (including the clover/SWterm). There 
are thus three independent improvement coeffi- 
cients, b q = 1 + 0(a s ), c' q — —0.25 + 0(a s ) and 
cngi = 0{a s ). The latter multiplies a gauge non- 
covariant term, but is nevertheless consistent with 
BRST symmetry. It turns out to vanish at tree 
level, but its one-loop value is not known — the 
one-loop calculation of the propagator in Ref. j| 
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only determines two linear combinations of the 
three coefficients. 

While, in principle, the three improvement co- 
efficients can be determine separately, this is diffi- 
cult in practice j5j. To understand this, note that 
the improved quark propagator is given by H 



Z q S = SL + ^ac' q + 2ac^GiSLii>—ab q mSL, 



(2) 



where Sl is the bare lattice propagator. To de- 
termine cngi and c' q one enforces, for large p, 



T*£(p) = 0, 



(3) 



since, in the continuum, chiral symmetry ensures 
that the r.h.s. of (|J) is asymptotically propor- 
tional to m/p 2 (up to logarithms). As can be 
seen from Eq. the c' q term gives a constant 
(p-independent) contribution to Tr S, which is to 
be chosen to cancel a corresponding contribution 
from Sl- To evaluate the cngi term, we can 
replace Sl by Z q S (since 0(a 2 ) terms are not 
controlled), and use the fact that, for large p, 
S = l/[i^£i(p 2 )]. Thus the cngi term is pro- 
portional to a/T, 1 {p 2 ). This is not a constant, 
and thus differs, in principle, from the c' q term. 
The two contributions are difficult to separate in 
practice, however, since Si(p 2 ) only varies loga- 
rithmically with p. 

Were one able to determine c' q and cngi, the re- 
maining improvement coefficient b q and the nor- 
malization Z q could be determined by enforcing 



Tr 7„S(p)~ 



i Tr 7„ j> , 
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where y, is a renormalization scale chosen so that 
sub- leading terms proportional to powers of m/fjL 
and Aqcd/a* can be ignored.^ If we decompose 
the inverse propagator as 



Sip)- 1 = i ^E 1 (p 2 ) + S 2 (p 2 )+0(a 2 ), 



(5) 



then Si(p 2 ) = 1 + 0{m 2 /fi 2 ). Note the pres- 
ence of sub-leading terms in Si, which are phys- 
ical, and will play an important role below. A 
possible definition of the improved, renormalized 
quark mass is then tor = ^(fJ 2 )- With these 
definitions, I note that 



s(pT 



= ii>+m R + 0{a 1 ,i>m 2 /p 2 ) . (6) 



Now I return to the main question: What er- 
rors are induced if we set cngi = 0, but keep all 
other improvement terms? Then we have an ap- 
proximately improved field 

-1/2, 
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I will denote the differences between true and ap- 
proximate coefficients as follows: 



Ab q = b q -b qi 



Ac' = c'-c 



etc. 



(8) 



The idea is then to adjust the approximate im- 
provement constants so as to compensate, to the 
extent possible, for the absence of the cngi term. 
In particular, I propose enforcing 



TrS(p) 



= 0, 



p 2 —(_L 2 ,m— 

which is a weaker version of Eq. (0), and 



= i Tr 7„ f , 



(9) 
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which has the same form as Eq. (Q). Since both 
conditions are satisfied by S, they apply as well 
to AS* = S - S. 

It is straightforward to work out the implica- 
tions of these new conditions using the result 

3 For later algebraic convenience, I have used a different 
convention for than the standard choice of Ref. Jlj], 
which involves derivatives. 



-[AlnZ°+mA6,+CNGi#]S , (11 



with cngi being the correct value of this improve- 
ment coefficient. Enforcing Eqs. (^,|l^) leads to^ 

Ac' q = Z°cxGi, AlnZ° = 0, Ab q = 2c NG im R /m . 

Since tor/to = Z m (l + ab m m) is of O(l) in per- 
turbation theory, as is Z q , we see that the errors 
in Cg and b q are proportional to cngIj and thus of 
0(a s ). This turns out to be generic for "unphys- 
ical" improvement coefficients, i.e. those that do 
not enter into on-shell matrix elements of gauge 
invariant operators. By contrast, there is no er- 
ror in Z q : ignoring cngi leads to the correct 0(a) 
improved normalization in the chiral limit. This 
is also generic: bilinears are 0(a) improved in the 
chiral limit even if cngi is ignored. 

Using the results for the errors in improvement 
coefficients, the error in the inverse propagator is 



A[S , - 1 ] = 2ac N GiS'" 1 (ii>+m R -S- 1 



(12) 



From this I conclude that, at generic large p, ig- 
noring Cngi leads to an O(a) error in the prop- 
agator. This is as expected since we are per- 
forming an approximate improvement. However, 
for p 2 — /i 2 , where S~ x = ij> + tor, it appears 
that AS^ 1 = 0(a 2 ). This would be very sur- 
prising, since it would allow the extraction of 
the exact improved renormalized quark mass — 
a physical quantity — in an approximate improve- 
ment scheme. 

The resolution is that one must include sub- 
leading terms in S _1 when evaluating Eq. (ft2|): 



Sip)- 1 



ii>[l + Oim 2 /p 2 )]+m R . (13) 



The subleading terms are physical, i.e. not lattice 
artifacts, but in perturbation theory they are of 
0(a 2 ) in Landau gauge Q. One then finds 

AiSip)- 1 ] cx ac N Gim 2 0(a 2 ) (14) 

which leads to 

Am R cx acNGiw 2 0(a 2 ) <^> Ab m cx CNGiO(a 2 ) . 

Thus the physical improvement coefficient b m 
does have an error if one ignores cngIi but this 



The first two equations were also obtained in Ref. 
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is of 0(a?s) and presumably much smaller than 
those noted above in unphysical coefficients. 

The only good news concerns £i(p 2 ). Since 
this is independent of m at large p, it can be eval- 
uated from the (correctly improved) propagator 
in the chiral limit. But in this limit one finds 
(from Eq. @) 

A[S(p,m = 0)- 1 ] = 2ac NGI p 2 [S 1 (p 2 ) - l]Si(;p 2 ) . 

Thus there is no O(a) error in Si, the j> part of 
S 1-1 , if it is evaluated in the chiral limit: 

£i(p 2 ) = Si(p 2 ) + 0(am) + 0(a 2 ) . (15) 

It is straightforward to extend the analysis to 
bilinears. The off-shell 0(a) improvement of bi- 
lincars involves only gauge invariant operators, so 
there are no other improvement coefficients analo- 
gous to cngi- In particular, one can work through 
the procedure for improving bilinears using am- 
putated correlation functions laid out in Ref . , 
but using the approximately improved propaga- 
tor S instead of S. I do not have space to present 
the details here, and report only the conclusions. 
As above, one must be careful to keep sub-leading 
contributions to physical vertices, or one can end 
up fooling oneself that 0(a) errors can be avoided 
at special kinematic points. 

For the scalar bilinear, off-shell improvement 
requires using 

S(x) = Z° s (l + ab s m)[qq(x) + ac' s E s (x)] , (16) 

where Es is an operator, defined in Ref. ||, which 
vanishes by the lattice equations of motion, and 
thus leads only to contact terms. On-shell im- 
provement requires knowledge only of bs, and the 
normalization constant Z s . I find that the error 
in these constants made by ignoring cngi is 

{Ac' q + Ac' s ) Z° s = 0x0!, AlnZ°=0, 

and A&5 0cc NG iO(a|) . (17) 

Combining these with the previous result for Ac' q , 
I find the same heirarchy of errors as deduced 
from the propagator analysis. The error in the 
off-shell improvement coefficient c' s is "large" , of 
0(a s ), while that in the on-shell coefficient bs is 
"small", of 0(a^), and that in the normalization 



vanishes. That Abs is proportional to a s and not 
a 2 is due to the fact that the sub-leading mft/p 2 
corrections to the scalar vertex vanish at one-loop 
order in Landau gauge. 

For all the other bilinears, the analysis is more 
complicated, because the determination of all the 
improvement coefficients requires applying one 
improvement condition at non-forward momenta. 
This introduces the dependence of vertices on 
momentum transfer, and, it turns out, leads to 
the generic expectation that Abr ~ Acr oc 
CNGiO(a s ), one power of a s larger than above. 
One also finds that A\nZ^ = 0. 

In summary, ignoring cngi leads to errors in all 
improvement coefficients, as expected, but these 
errors are quite small for those needed for calcu- 
lations of physical matrix elements. It may be 
that other sources of systematic error, such as 
the need to subtract large 0{a 2 p 2 ) corrections, 
exceed the error made by ignoring cngi- Further- 
more, leaving out the cngi term has no effect on 
the normalization constants of the operators in 
the chiral limit. 

On the other hand, it is preferable to have a 
method without uncontrolled errors. How can 
this be acheived using NPR, given the difficulty in 
determining cngi? One method is to input knowl- 
edge of one or more improvement coefficients ob- 
tained by other methods, e.g. that based on Ward 
identities. For example, if one knew b m , one could 
adjust cngi until the result for wir had the cor- 
rect ratio of quadratic to linear dependence on the 
quark mass m. This seems to be equivalent to the 
method suggested here by Bhattacharya ||. 
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